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Abstract 

If collider experiments demonstrate that the Minimal Supersymmetric Standard Model 
(MSSM) is a good description of nature at the weak scale, the experimental priority will 
be the precise determination of superpartner masses. These masses are governed by the 
weak scale values of the soft supersymmetry (SUSY)-breaking parameters, which are in 
turn highly dependent on the SUSY-breaking scheme present at high scales. It is there- 
fore of great interest to find patterns in the soft parameters that can distinguish different 
high scale SUSY-breaking structures, identify the scale at which the breaking is commu- 
nicated to the visible sector, and determine the soft breaking parameters at that scale. 
In this work, we demonstrate that 1-loop Renormahzation Group (RG) invariant quanti- 
ties present in the MSSM may be used to answer each of these questions. We apply our 
method first to generic fiavor-blind models of SUSY-breaking, and then examine in detail 
the subset of these models described by General Gauge Mediation and the constrained 
MSSM with non-universal Higgs masses. As RG invariance generally does not hold beyond 
leading-log order, wc investigate the magnitude and direction of the 2-loop corrections. 
We find that with superpartners at the TeV scale, these 2-loop effects are either negligible, 
or they are of the order of optimistic experimental uncertainties and have definite signs, 
which allows them to be easily accounted for in the overall uncertainty. 



Introduction 



The Standard Model (SM) provides an excellent description of the experimental data collected 
at high energy physics facilities. However, it has a series of shortcomings and is believed to 
provide only an effective theory at energies below a scale of the order of 1 TeV. The Minimal 
Supersymmetric Standard Model (MSSM) addresses most of the SM shortcomings: it provides a 
solution to the hierarchy problem, is consistent with the unification of couplings at high energies, 
leads to an understanding of the origin of the negative Higgs mass parameters generating the 
breakdown of the electroweak symmetry, and includes a natural Dark Matter candidate [IH3]. 
Furthermore, the MSSM provides a weakly interacting extension of the SM in which the SM-like 
Higgs mass is bounded from above by about 130 GeV and will be within observational 

reach of colliders within a few years p^HSOj . Finally, MSSM sparticle effects on the precision 
electroweak observables efficiently decouple for sparticle masses above a few hundred GeV, in 
agreement with the best fit to these measurements (see, for example, [5TH55] ). 

If sparticles corresponding to the MSSM are discovered at the LHC or any future collider, 
resolving the properties of the supersymmetry (SUSY)-breaking mechanism will become an 
important theoretical and experimental problem. Since most well-motivated scenarios place 
the "messenger scale" at which SUSY-breaking is communicated to the MSSM far beyond 
direct experimental access, it is likely that the renormalization group (RG) will be our primary 
tool for investigating the nature of this mechanism. 

The standard approach is "top-down," where a high scale model with a limited number 
of parameters (for instance the CMSSM or mSUGRA) is assumed, and the parameters are fit 
to the low scale data using RG evolution. This method is well studied, both in the context 
of present constraints from dark matter searches and low-energy observables as well as for 
future LHC projections (for a sample of references, see [MlflT] ). It is now quite a sophisticated 
technique (see for instance ^j), but it has a few limitations. Primarily, the reliability of the fit 
may be reduced as the number of high scale parameters is increased. It is also subject to the 
behavior of uncertainties under RG flow: if they tend to spread as the scale increases, the fit 
will become weaker. The top-down approach is also sensitive to the experimental uncertainties 
in gauge couplings and SM particle masses, and multiple iterations of RG evolution from 
high to low scales are required to consistently incorporate SUSY corrections to the Yukawa 
couphngs |^3] - H6] . 

Alternatively, a direct "bottom-up" approach has been studied [iTHST] and found to pos- 
sess utility complementary to the top-down method. The bottom-up reconstruction has the 
advantage of being easy and quite transparent: simple analytic formulae are used to convert 
measured pole masses into low scale values of the Lagrangian parameters, including dominant 
radiative effects. The running parameters are then evolved up to the high scale at which the 
mediation of SUSY-breaking takes place and their structure can be analyzed. This method is 
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relatively fast and attractive, and it is less sensitive to the total number of free parameters, 
but it also presents challenges: foremost, one does not know the scale that should be used 
as the ultraviolet boundary. Additionally, uncertainties in the running parameters must again 
be propagated carefully under RG evolution. As in the top-down approach, the /5-functions 
for each sfermion mass are subject to the present experimental uncertainties in the gauge and 
Yukawa couplings. Furthermore, these /3-function depend on many MSSM parameters even 
at 1-loop order, coupling all of the soft scalar masses to each other through the hypercharge 
D-term, as well as to the gaugino masses and the soft trilinear couplings. In Ref. [52|, it was 
demonstrated that the failure to infer even one soft mass experimentally may be sufficient to 
drastically alter the high-scale prediction for other masses. 

Here we would like to suggest an alternative technique for the bottom-up reconstruction 
program. Within the parameters of the MSSM, neglecting first and second generation Yukawa 
couplings, there exists a significant set of combinations of the parameters that are RG invari- 
ant (RGI) at 1-loop. We propose that these RGIs may be used to efficiently test and extract 
parameters of high scale SUSY-breaking models, including in certain cases the messenger scale. 

There are several reasons why RGIs are of particular interest. The first is that in some appro- 
priate bases, subsets of the RGIs are predicted to vanish in different classes of SUSY-breaking 
models. A low scale measurement of any of these invariants may therefore provide immediate 
evidence for or against several attractive high scale models without requiring knowledge of the 
value of the high scale. The second reason is that for broad classes of SUSY-breaking models, 
most or all of the fundamental high scale parameters can be expressed uniquely in terms of 
the nonzero RGIs. If measurements of the first type of RGIs indicate that a particular class 
of models is favored, the second type may be used to constrain the parameter space of those 
models. Thirdly, none of the RGIs individually involve all of the parameters of the theory. 
This property limits the set of physical masses and couplings that must be measured before 
high-scale theories can be tested. Finally, from a practical point of view, the use of RGIs 
simplifies the standard bottom-up approach: it avoids the integration of RG equations and the 
complicated propagation of errors between scales. 

Although RG invariance only holds strictly at leading-log order, 2-loop corrections are ex- 
pected to be quite small and are likely to be negligible compared with experimental uncertain- 
ties. To confirm this expectation, we analyze the 2-loop effects and demonstrate that they can 
either be neglected, or absorbed to a good approximation into a simple shift of the measured 
value of the RGIs. For each RGI, we show that the shift can be computed from a linear function 
of a few dominant sparticle masses. 

In this work, we will study RGIs within theories satisfying the following conditions: 

• The effective theory at the electroweak scale is the MSSM; 

• No new physics alters the 1-loop MSSM /3-functions below the messenger scale, at which 
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SUSY-breaking is transmitted to the observable sectoiQ 

After deriving a set of invariants under these assumptions, we will apply our method to 
test the class of theories in which the SUSY-breaking mechanism is flavor-blind. In addition 
to the general case, we will consider two major subclasses of flavor-blind theories. The first 
is general gauge mediation (GGM), where the interactions between the SUSY-breaking sector 
and the sfermion/gaugino fields of the MSSM are controlled entirely by the MSSM gauge 
couplings [M]- Integrating out the hidden sector generates soft mass terms for these fields that 
can be parameterized by six constants. Soft masses are also generated for the Higgs doublets, 
but it is likely that additional, non-gauge couplings link the Higgs and hidden sectors. We 
demonstrate that certain RGIs may be used to test the GGM hypothesis, and if the data is found 
to be consistent, other RGIs may be used to extract the full set of GGM parameters, including 
the soft Higgs masses at the messenger scale. The second type of SUSY-breaking structure we 
examine is the Constrained MSSM (CMSSM) with non-universal Higgs masses (NUHM). In this 
case, fiavor-blindness is maintained with fewer input parameters than in GGM. Consequently, 
all messenger scale parameters may be expressed in terms of a subset of the RGIs. Among the 
remaining invariants, several nontrivial relationships must be satisfied for the low energy mass 
spectrum to be consistent with the CMSSM + NUHM. 

The RGI quantities discussed in this work have all appeared previously in the literature in 
various guises, but either as limited subsets or not in the context of a bottom-up determination 
for the high scale SUSY-breaking parameters. In Ref. [52] several invariants were introduced for 
the purpose of testing high-scale fiavor universality in the presence of extra intermediate-scale 
GUT multiplets. Ref. [22] also commented on the use of two quantities, Dy and Db^l (defined 
in Section 2), for probing GUT embeddings of the MSSM. These functions are related to the 
RGIs but are not themselves invariant in general: we will discuss the use of the corresponding 
RGIs in greater detail for testing gauge mediation and the CMSSM. In Ref. [HS] most of 
the invariants we will use are listed, but in a basis without the properties discussed above, 
rendering them less relevant for probing high scale models with low, single-scale measurements. 
In Ref. [SI], Dyi Db-Li and one other invariant combination of them were listed as potential 
tests of GGM; they will be included and expanded upon in our discussion. Finally, Ref. [561460] 
discuss in detail some of the invariants in the context of sum rules for mSUGRA and a variety 
of SUSY-GUTs. 

We attempt to go beyond previous work on the subject of RGIs by providing a full set of 
invariants in a useful basis, a program for using them to explore the structure of the messenger 
scale effective theory, a discussion of experimental uncertainties, and an analysis of 2-loop 
effects. Our presentation is divided into four sections. In the first, we will derive the RGIs 
in a convenient basis from symmetry arguments. In the second we will discuss tests of flavor- 

^Here we neglect the possibility that strong couplings in the hidden sector could affect the running of the 
soft scalar masses, as pointed out in |53j . 
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blindness, and under that hypothesis we will express the high scale soft parameters as simple 
functions of RGIs and a single gauge coupling at that scale. We will then use the RGIs to test 
and determine the parameter spaces of GGM and of the CMSSM + NUHM. In the third section, 
we analyze experimental prospects and 2-loop corrections for the RGI method in the context 
of GGM. We discuss briefly the minimal set of physical observables necessary to determine the 
soft parameters in the RGIs, and estimate the precision with which the RGIs may be used 
to constrain the high scale parameter space given a set of uncertainties in the low scale soft 
parameters. We derive simple estimates that may be used to subtract out the bulk of the low 
scale 2-loop shifts in the RGIs, and demonstrate that residual 2-loop contributions are typically 
small compared with optimistic assumptions for experimental uncertainties in the soft masses. 
We reserve the final section for our summary and conclusions. An appendix collects some 
analytic approximations for the 2-loop /3-functions. 



1 1-loop RGI Combinations of Soft Masses and Gauge Couplings 

In this section we will derive fourteen RG invariant^ and relate several of them to symmetries 
of the MSSM Lagrangian. For reference, the RGIs are collected in Table |2] 

We begin with RGIs constructed solely out of soft scalar masses. We will assume that the 
soft sfermion masses are flavor diagonal in the super-weak basis and that the first and second 
generation masses are degenerate at the input scale, as is strongly motivated by low-energy 
flavor mixing constraints. Similarly, we assume that there are no new sources of CP-violation 
in the soft SUSY-breaking sector beyond those induced by the Yukawa couplings. We will also 
neglect first and second generation Yukawa and trilinear couplings since they give very small 
contributions to the evolution of the soft SUSY-breaking parameters, smaller than the 2-loop 
corrections associated with the gauge and third generation Yukawa couplings. In general, these 
approximations tend to be excellent for models in which there exists flavor universality at the 
messenger scale, and all flavor violation effects are due to radiative corrections induced by 
Yukawa couplings. Finally, we assume that the right-handed neutrino, if it exists, effectively 
decouples from the spectrum and the RGBs (see discussion in Section 3). With these conditions, 
the soft scalar masses rrii obey the following RGEs at 1-loop: 



IQvr^^ = E yW(^' + ^? + + ^^k^"'') - 8 J2 Ca{^)glm' + hglDy, (1.1) 



dt 

jk 



few more RGIs exist, such as I2 and I4 of Ref. 55 . However, they are not needed in our investigation of 
flavor-bhnd theories, and they require more MSSM parameters than the set considered in this work, including 
the soft trihnear couphngs and ^. 
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Table 1: U(l) Representations in the MSSM 



Particle 


Y 


B 


L 






1/6 


1/3 







V ' J L 


-1/2 





1 


u 




-2/3 


1 /in 

-1/3 





d 




1/3 


-1/3 





e 




1 





-1 




Hu 


1/2 












-1/2 









where 



Dy = Tr(rm^) 

= 5Z ("^^~2m? +mj-m| + m?) +m^^ -m^^. (1.2) 



gen 



The first sum in Eq. (1.1 ) is over all the degrees of freedom available to run in self-energy loops, 
the second sum is over the three gauge groups a of the MSSM, Ma and Aij^ are the gaugino 
masses and soft trilinear parameters, Ca{i) is the quadratic Casimir for the representation i of 
a, and t = \og{^/Mz). The trace in the definition of the hypercharge D-term, Dy, runs over all 
chiral multiplets. A summary of the B, L and the hypercharge representations for the particle 
content in the MSSM is given in Table [TJ 

We first construct linear combinations, Di, of the soft masses that evolve only with Dy, 

A = Tr(g,m2), ^ = aiDy, (1.3) 

where Qi and are constants, and again the trace is defined over all chiral multiplets. With 
five soft SUSY-breaking masses contributed per generation (for a total of ten soft masses if 
the first generation masses are identified with those of the second generation) and the two soft 
SUSY-breaking Higgs mass terms, twelve coefficients Qi are associated with each combination 



Di. That there are six independent satisfying Eq. (1.3) can be seen from Eq. (1.1 ) as follows 
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For the Yukawa terms to vanish from the /3-function, the Qi must correspond to charges of a 
global symmetry of the classical Yukawa potential. This affects only the third generation in our 
approximation, so it implies three independent constraints on the twelve Qi. For the gaugino 
terms to cancel, the symmetry must have vanishing mixed anomalies with the SM gauge groups. 
This supplies three more independent constraints on the Qi. 

Furthermore, we can construct a basis in which five of these six combinations will also satisfy 
TiQiY = 0. This condition eliminates the Dy term from the RGEs. These five combinations 
are then genuine RG invariants, independent of the vanishing of Z)y. 

Let us first consider baryon number, Q = B, and lepton number, Q = L. Classically, these 
are symmetries of the action, but both are anomalous in the MSSM. Within our approximation, 
both the baryon and the lepton number associated with each generation is separately conserved 
at the classical level and their anomalies are flavor-independent. Hence the differences between 
the first (or second) and third generation baryon and lepton number are anomaly-free. With 
this in mind, one can define RGIs associated with the new global symmetries: Bi^ and L13 j52], 

Dbi3 = ^Bi - Db-3 

= 2m|^ - ml - ml - 1m\^ + ml^ + m\ , 



= 2m\-m\-2m\^^m\, (1.4) 

where numbers in the subscripts indicate the generation. 

We can also consider non-anomalous f/(l) symmetries. The obvious choices are hypercharge 
and [B — L). However, since Try^ and TtY{B — L) are non- vanishing, these combinations will 
evolve with Dy. 

In the hypercharge combination, the Dy dependence can be eliminated by using nearly the 
same trick with family non-universality as we did with B and L. In this case, we must include 
the Higgs doublets with the third generation, since their soft mass evolution is linked by the 
Yukawa couplings. The proper RG invariant combination is given by 



10 
13 



For {B — L), generation subtraction is redundant, since we can construct it out of Db^^ and 
Di^g. However, even restricting the trace to a single generation, Dy and Db^l evolve only 
with Dy. Therefore we can construct a fourth RGI combination which depends only on the 
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soft masses of the first generation j54j , 



= —6m% — 3m? + 9mj + Gm? — m? . (1.6) 

(q/l "1 ai 1^1 ^1 ^ ^ 

Here xi is used to indicate that this charge assignment is related by an overall re-scaling to the 
U{1)^ symmetry generated by the spontaneous breaking of Eq. 

The Eq breaking generates an additional U{1), Eq — )• SU{5) x U{1)^ x U{1)^, but the 
[/(l)^ symmetry is anomalous when restricted to the MSSM and involves the Higgs sector, 
which prevents a family non-universal invariant combination. Fortunately, by taking the linear 
combination (the "inert f/(l)"' [SI]), 

^um, + fu(i),. (1.7) 

and setting the charge of the first generation left-handed sleptons to zero, we obtain an anomaly- 
free symmetry, U{l)z, and the corresponding RGI |52] : 

Dz ^ 3m| + 2m|^ - 2m^^ - 3m| . (1.8) 

This new symmetry is evidently independent from those already discussed, as it depends only 
on 171% and not on the combination m?j — m?T . 

Of course, the sixth combination must now run with Dy- The RGE for Dy is given by 

16vr — = -g,Dy, (1.9) 

so this combination can be taken to be Dy itself. Db~l also evolves only with Dy, and as 
mentioned previously, so do the restrictions of Dy and Db-l to any given generation (plus 
Higgs in the case of the third generation Dy). These quantities are then RGIs if Dy vanishes, 
although only those related to Dy are linearly independent of the invariants introduced above. 
However, one can construct a genuine independent invariant out of Dy that is nonlinear in 
running parameters. Recall that the gauge couplings obey homogeneous RGEs at one loop: 

levr^^ = g'r (Tr„/,(n) - SCriG)) , (1.10) 

where Ir{n) is the Dynkin index of the chiral multiplet n, and Cr{G) is the quadratic Casimir 
invariant of the adjoint representation. The gi group theory factor is equal to 33/5 in the 
MSSM, and so we obtain the RGI, 

lYa = ^. (1.11) 

9i 
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Table 2: 1-Loop RG Invariants in the MSSM 



Invariant 


Symmetry 


Dependence on Soft Masses 




Bi — B3 


2(m| - m| ) - ml + ml - m| + m| 




Li — L3 






Xi 


3 (3m J — 2(m^ — m| ) — m? ) — m? 








Dz 


z 




Iva 


Y 


(mjj^ - mjj + J2geniml - '2ml + - m| + m?)) /gl 


















M| + ^ (9(m2 - m|J + 16m| - m?J 












1/^? - 33/(5^i) 


^93 




+ 33/(15^1) 
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Additionally, because of the form of Eq. (1.10), we obtain two additional RGIs 

33 



•■92 



1 33 



(1.12) 



whose values may be computed from measurements at the scale Mz- In the MS scheme, for 
instance, they are approximately —10.9 and 6.2, respectively. 

We now construct RGIs with gaugino mass dependence. The three soft gaugino masses, 
Mr, evolve with 1-loop RGEs given by 



IQ-n'^dtMr = glMr (2Tr„/,(n) - QCr{G)) . 



:i.l3) 



Again replacing in Eq. (1.13) with dt\og{gr) and the group theory factor, we recover the 
well-known set of three RGIs, 

Jb,. = Mr/gl . (1.14) 

Finally, there are three RGIs that mix gaugino and sfermion soft mass parameters. Following 
the method of Ref. |55], these can be obtained easily from the 1st generation masses because 
of the absence of Yukawa couplings and the homogeneity of the and Dy RGEs, 



(1.15) 

The invariant Imi involves the subtraction of potentially large squark masses squared with 
a large coefficient. In models possessing a large mass hierarchy between strongly and weakly 
interacting particles, the splitting between m^^ and ttlu^ is likely to be much smaller than either 
mass, whereas the experimental uncertainty in Im-^ grows approximately linearly with both mj^ 
and rrtu^ . Therefore, the measured value of this RGI is likely to be consistent with zero within 
experimental errors. A similar effect may occur with i but in this case the squarks coefficient 
is significantly smaller. One can form a combination of Im^ and that is independent of the 
squark masses, and therefore suffers lesser uncertainties: 

= Ml+^-^{?>Ml + 2ml^+mlY (1.16) 

In most cases, /mi2 may be determined much more precisely than Im^ or Im2- Indeed, since all of 
the terms are positive, the percentage error in its determination is governed by the percentage 



Im2 



= Ml 



33 



= M^ + 



Mi 



24 
3 

16 



2 2 2 

di "1 ei 



9{mr -m-J + 16m. - 



5m| + m?^ 



m 



ei 
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error in the measurement of the weakly interacting sparticle masses. 

With the exception of lyai the RGIs given in Table 1 of Ref. [531 may be obtainec^as linear 
combinations of those we have listed in Table [2j However, the explicit relationship to symmetries 
in our basis allows a clean application of the invariants to the problem of reconstructing the 
theory at the messenger scale. 

2 Extraction of High Scale Parameters from RGIs 

We now examine the extraction of high scale SUSY-breaking parameters from the RGIs con- 
structed in the previous section. The five soft masses contributed by the degenerate first and 
second sfermion generations, five soft masses from the third sfermion generation, two Higgs soft 
mass parameters, three gaugino masses, and the gauge couphngs (assuming no knowledge of 
the high scale) constitute eighteen unknowns (more are introduced if the right-handed neutrino 
or other new physics is included). In the most generic case, assuming all RGIs are non-zero at 
the high scale, fourteen of these parameters may be obtained. Although the system is under- 
constrained, the high scale degrees of freedom are considerably reduced by the determination 
of the invariants. 

As mentioned previously, a large class of SUSY-breaking models favor flavor-universal in- 
teractions. Even without the ability to extract all high scale parameters, using the RGIs to 
test flavor blindness would be a non-trivial result. Furthermore, in the most generic flavor- 
blind models, the number of unknown parameters is reduced to thirteen, whereas the number 
of nonzero RGIs is twelve. We will consider the generic case below, as well as the specific 
models of GGM and the CMSSM with non-universal Higgs masses. The high scale values of 
the invariants listed in Table [2] are tabulated in Table |3] for these three models. 

2.1 Generic Flavor-Blind Models 

The most immediate consequence of flavor blindness is the vanishing of D and D^^^ . There- 
fore these invariants provide us with a direct test of the flavor-independent hypothesis with 
a minimal set of measurements. More precisely, they allow this hypothesis to be ruled out: 
measuring Db^^ 7^ or -Dlis 7^ at the low scale implies high-scale family non-universality; 
however, as noted in Ref. |52j, measuring Db^^ =0 and Dl^.^ = at the low scale does not 
necessarily indicate high-scale universality. 

Current experimental data from flavor physics strongly motivates a flavor-universal me- 
diation mechanism for SUSY-breaking (for a selection of Minimal Flavor Violation studies, 
see [621479] ). Accordingly, if Db^^ and Dl^^^ are found to vanish, it is reasonable to proceed a 

few typos appear in Ref. [55]: the coefficient of A/| in I12 should be +3/2, and similarly for /17. The 
coefficient of should be 1/33, and all coefficients of S should be multiplied by 13/33. Note that with these 
corrections, five invariants in the list /5-/17 of Ref. |55j can be expressed as linear combinations of the others. 
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Table 3: Values of the 1-loop RGIs in terms of high scale soft parameters. 



Invariant 


Generic Flavor Blind Model 


GGM 


CMSSM with NUHM 


























3(3m|-2(m|-m|)-m?)-m? 





5ml 




^ f 3(m| - 2m| + m| - m| + m|) 
+ \^{m\ - mjjj) 






Dz 


2(m|-my 








^3(m^ — 2m| + m| — m| + m|) 

- + m%J hi 


{K - Sd) /gl 


{K - Sd) /gl 






MBr 


m,/2/gl 




^I' + f (^i-^l+^i) 


gf {{MB,f + |A) 


^1/2 + f ^0 






gt {{MB^f + \A^) 


^1/2 + Wo 






g\ {{MBsY - lA,) 






w -10.9 


w -10.9 


a; -10.9 


^93 


p=i 6.2 


^ 6.2 


f=i 6.2 
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step further and attempt to extract constraints on the high-scale values of flavor-blind MSSM 
soft parameters from the RGIs. 

The ten free soft mass parameters can be expressed uniquely in terms of the ten invariants 
D^-^ through Im^ listed in Table |2] and are given in Eq. (2.17)-(2.18). Note that these relations 
depend on the three gauge couplings at the high scale; with the invariants Ig^ and Ig^ defined 
in Eq. (1.12), they may be expressed entirely in terms of gi. Equivalently, one can reduce the 
degrees of freedom at the high scale to a single parameter which can be taken to be the value 
of that scale. In particular this permits tests of more restrictive flavor-universal models such as 
mSUGRA, taking gi at the GUT scale. However, if the high scale is not known, one may take 
advantage of the fact that the gauge couplings do not vary too wildly with scale and simply 
guess values for them, at the cost of introducing a further uncertainty. We will discuss this 
possibility further below. 



Ml = gllB, , 
M2 = gllB, , 

M^ = gllB,, (2.17) 



2 2 ^ -'■3 n , 9l T 

ml = ^(26l^n3H + 11^X1 - H^aUk - hh) - gliva)) , 

ml = -^{78DYr,, + 33D^, + 20(4(((?^4 - Uglily) - {Im, - IUm,)) + ^gflya)) , 

mj = Y^i^y,,^ + 33/^x1 - H^iia'Jk - ^^atilj - [hi, - 44/M3)) - Hiya)) , 

m\ = ^— (78Dy,,„ - 627L>., 
3960^ 

- 20((^t4 + 297^72^4 - ^7^9pk) - {hi, + 297/m, - 176/^3) - ^gflya)) • (2.18) 
In the above, all couplings and soft parameters are assumed to be given at the messenger scale. 

2.2 General Gauge Mediation 

Gauge mediated SUSY-breaking encompasses a broad class of models in which flavor blindness 
is perhaps most natural (SDHHS]- R-^f- [51] General Gauge Mediation (GGM) was defined 
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as any theory in which all SUSY-breaking effects decouple from the MSSM in the limit of 
vanishing MSSM gauge couplings. 

Let us first review the formulae for the soft masses in GGM at the messenger scale given in 
Ref. |5l]. The soft sfermion masses can be parameterized by a set of three parameters, ^rj^ 



3 



m% = J29tCr{f)Ar, (2.19) 



where the sum runs over the gauge groups of the MSSM, and the Ar can be expressed as loop 
integrals over hidden sector current-current correlation functions, 

(2.20) 

The gaugino masses are expressed in terms of three more parameters, MBr, given by 

Mr = glMBr , (2.21) 

where M is the scale of SUSY-breaking and 

MBr = -^e"^(j„(0)j^(0)) . (2.22) 

In order to generate a Higgsino mass parameter, fi, and soft term, 5^, of the correct order, 
gauge mediation may need to be supplemented by additional SUSY-breaking contributions in 
the Higgs sector. Therefore, we assume that in the case of the soft Higgs masses, the expression 



given in Eq. (2.19) may be modified. 



<,=<,+^<i- (2-23) 



With these definitions for the soft masses, we can now calculate the high scale values of the 
RGIs in GGM, listed in Table [3j By inspection, we see that they can be separated into 

• Three that test consistency of GGM; 

• Six that probe the high-scale mass parameters of pure GGM; 



We neglect a possible explicit hypercharge Fayet-Iliopoulos (FI) dependent term, gfYfS^, in Eq. (2.19). The 



FI term may drive some of the rrij to negative values, and we assmiie it is fixed to zero by a Z2 symmetry in the 
hidden sector. For further discussion of difficulties associated with FI terms, see [521 [HZ]- We comment again 



on this term below Eq. (2.28) 
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• Five that are sensitive to extra structures linking the Higgs and SUSY-breaking sectors 
and to the values of the gauge couplings at the messenger scale. 

The flavor-independence of gauge mediation is manifest in the formulae listed for the soft 
masses, and therefore, as before Db^^ and D^^^ should vanish at all scales. 

The third invariant testing consistency of data with GGM is D^-^ [54j. It should also remain 
approximately zero at all scales, which can be seen as follows. The condition for anomaly 
cancelation between a U{1) symmetry Q and another gauge group G with generators in 
the representation i reads 

Trgt5t5 = 0, (2.24) 

where the trace runs over all degrees of freedom. If Q commutes with all generators of G, this 
can be written equivalently as 

= 5^Q,Trt^(z)t5(z) , 

i 
i 

= TiQ,CG(t), (2.25) 

where the sums run over all representations i of G and d{i) is the dimension of i. The sfermion 
masses in GGM depend on the gauge couplings and the Casimirs of the gauge groups; hence, 
if the mixed anomalies between Q and the MSSM gauge groups cancel, we can multiply both 



sides of Eq. (2.19) by the charges Q and take the trace to obtain 

Dq = TiiQm^) = J2 9%Ag Tr Q,Cg{i) = . (2.26) 



The second equality in Eq. (2.26 ) is only valid if the Higgs fields are neutral under Q; otherwise, 
the modifications of Eq. (2.23) imply Dq = 2{Qh^5u + QH^^d)- Choosing Q = xi, the Higgs 
fields are uncharged, and thus vanishes. This property may differentiate GGM from alter- 
nate theories with universal soft masses such as mSUGRA, where does not vanish unless 
a relevant right-handed neutrino is included. 

This argument also suggests Dy and Db-l vanish and are RG invariant, since both are 
anomaly-free symmetries of the MSSM and both run with Dy- However, as mentioned above 



the modifications to the Higgs sector in Eq. (2.23) may spoil the vanishing of Dy. Furthermore 



although Dy is zero in pure GGM, it also vanishes in other SUSY-breaking scenarios, including 
mSUGRA. 

These first three RGIs are sufficient to either rule out GGM or demonstrate consistency 
of GGM with the data. However, they do not constrain the (A^, MBr, Su, Sd, M) parameter 
space. Therefore, we now turn to the RGIs with non-zero values at the messenger scale. 
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In the most generic case where Dy (and thus the invariants Dy^zh ^^"^ Iyo) is nonzero, we 
can determine the gauge couphngs at the high scale. The following relations are satisfied at 
the input scale: 



Dy: 



10 



l-Ya 



13 

9\ 



(2.27) 
(2.28) 



Therefore, if we can calculate the value of Dy^^hI ^Ya at the low scale, we can infer the value of 
g\ at the input scale via 



10 / 



Ya 



(2.29) 



Then the Ig^ RGIs given in Eq. (1.12) can be used with Eq. (2.29) to compute g2 and at the 
high scale. 

The invariants Dz and Dy^j^^ can be used to extract the individual input-scale values of 5u 
and 5^: 



2 V r 

2 



(2.30) 



Now we turn to the RGIs with explicit dependence on the gaugino mass parameters to 
extract information about the Ar and MBr- From the Ig^ we immediately obtain 



MBr 



I Br 



(2.31) 



For the sfermions, using the /m^, we can obtain the Ar 



10 fhl, r2 



A,. 



A, 



2 fh'ii j2 



3 V gi 



(2.32) 



where the g^ are the gauge couplings at the unknown scale, M, that can be deduced from 
Eq. (2.29) (note that the gr used to compute the I Br are at the low, measurement scale)[^ 
If 5u is found to be equal to 5^, and thus the gauge couplings at the high scale cannot be 



^If there is an explicit hypercharge FI term at the input scale, the only effect on the RGIs is to change the 
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extracted from the RGIs, explicit RG evolution may be the only available method for deter- 
mining M, choosing it to be the scale at which = "^Ig- This would require knowledge not 
only of all up-type soft mass parameters at the low scale, but also of the soft trilinear coupling 
At of the stop to the Higgs, which appears in the m|g RGE|^ Moreover, as discussed previously, 
the evolution of the soft masses is linked to the evolution of all other parameters via the Dy- 
terms. Therefore, it is expected that there will be large uncertainties in the determination of 
the messenger scale by this method. 

As mentioned before and as we will see in Section 3 below, Im^ and could suffer from 
large uncertainties induced by cancelation between potentially large masses. If Imi2 is used 
instead, this pitfall is avoided, but in general only the following correlation between Ai_2 and 
gi^2 can be obtained, 

33 1 1 

gfj^A, + gt-A, = Im,, - g\ll^ - llgtll, . (2.34) 

In the above, the coefficients of the Ib^ and Ai terms are significantly smaller than those of 7^2 
and yl2, particularly for low values of the messenger scale where g2 — 2gl. Therefore, unless 
Ai ^ A2, for low messengers scales, the above equation then gives the approximate expression 
for A2, 



M^^A^- ) . (2.35) 



2 

2.3 Constrained MSSM with Non-Universal Higgs Masses 



The CMSSM + NUHM [HSHM] is another common model realizing fiavor universality. The 
messenger scale is the GUT scale, however, here we will consider the more generic possibility 
of M 7^ Mgut, as for example occurs in Mirage Mediation [951497] . The sfermions are given a 
common soft mass mo and the gauginos share a soft mass mi/2 at the scale M. As in the GGM 
case, the Higgs masses are allowed to differ from the sfermion masses by the factors 5„ and b^- 
From the fourth column of Table [3} we see that the non-vanishing RGIs overconstrain the 
(mo, mi/2, gr{M), 6u, Sd) system, allowing multiple ways of testing consistency and extract- 
ing parameters. Both experimental errors and 2-loop contributions to the RGIs should be 
considered to decide which avenues to use for each purpose. The unification of gaugino mass 



input scale value of lya ■ 

Su^^Su_Sd^ (2.33) 



Eq. (2.32) for the A,, remains unchanged; however, this shift is sufficient to spoil the extraction of gi at the 



input scale from the inhomogcneitics in the Higgs sector. 



^One could equally seek the scale where m?- — m| 
be more difficult to determine experimentally than At- 



One could equally seek the scale where m - — m - , but it is likely that the soft trilinear coupling Ab will 
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parameters at some scale M* demands the following consistency relationships: 



T 33 



(2.36) 



33^ 



(2.37) 



where mi/2 is now a universal gaugino mass parameter at M*, which is not necessarily con- 
strained to be equal to M by these relations. 

In the MSSM the gauge couplings unify at the grand unification scale, and therefore at any 
arbitrary scale the following relationship is satisfied 



'93 



-1/ 

Y 92 ■ 



(2.38) 



Together Eqs. (2.36)-(2.38) indicate that gaugino mass unification requires the simple condition 



^/b, -/B, + ^/ij3 = 0. (2.39) 
This relation is necessary but not sufficient to ensure unification. In addition, a second require- 



ment following from Eqs. (2.36)-(2.38) is 



51, 



(2.40) 



92 



If we demand that Eq. (2.40) should be satisfied for a sensible value of M*, between about 10" 



and 10^^ GeV, then the Is. should obey the inequality 



Ibi ^ 5Ibi — 33/^2 ^ Ibi 



5/, 



92 



In the particular case where the gaugino masses unify at M* 



M. 



GUT, Jfii 



'52 



(2.41) 



Ibs and 



Eq. (2.39) holds directly, while the upper bound of Eq. (2.41) is approximately saturated. Note 



that it is not unreasonable for a model to obey Eq. (2.39) while still violating Eq. (2.41). 



This happens, for example, in anomaly mediated SUSY-breaking [98, 99J, where the Ib^ are 



proportional to the /3-function coefficients of the gi, and so Eq. (2.39) holds. The gaugino 



masses do not unify in anomaly mediation, however, spoiling Eq. (2.41) (and the extracted 



value of mi/2 from Eq. (2.36) or Eq. (2.37) would be zero). Note further that these relations 



involve invariants sensitive only to the gaugino sector, and therefore are expected to have small 
experimental uncertainties. 

Eqs. (2.36)-(2.41 ) guarantee gaugino mass unification, but do not enforce M = M*, which 



is a requirement in the CMSSM. In mirage mediation or minimal gauge mediation, the masses 
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unify, but M ^ M* in general. Therefore, to differentiate tlie CMSSM from these other models, 
one may test the equality of M with M*, if 5u 7^ by equating the messenger scale value of 
gl calculated via Eq. (2.29) with the M* value of gl found from Eq. (2.40). Additionally, the 
extracted value of mi/2 from Eq. (2.36) has to be consistent with the gaugino mass calculated 
from the other RGIs. These considerations enforce that the following set of relations, in addition 
to Eq. (2.39), have to be satisfied for the spectrum to be consistent with the CMSSM: 



2(5/bi -33Jb2^ 

- 33/B2 \ ^ 33L> 



5/, 



92 



+ 



Xi 



40 



57 



92 

- 33/B2 A " 3D 



5/, 



92 



Xi 



16 



(2.42) 



If the data satisfies Eq. (2.39) and Eq. (2.42), one can proceed to the extraction of parameters. 



D^^ immediately yields the value of the soft sfermion mass mo, and Eq. (2.30) can still be used 



to obtain 5^ and 5u- Turning to the 2-loop /3-functions for the RGIs listed in Appendix A, we 
find that /jv/i is a 2-loop invariant in the approximation that gi, the soft trilinear couplings, 
and the lepton Yukawa couplings are set to zero. It may therefore be useful to extract mi/2 
from Jjvfi. However, the advantage of 2-loop invariance can be mitigated by the consideration 
of experimental errors, which as noted previously may be large for Im^- In that event mi/2 



should be taken directly from either Eq. (2.36) or Eq. (2.37), for which the experimental errors 
are expected to be small: 



D 



rria 



XI 

5 ' 



mi/2 



5/bi — 33/^2 



5/, 



(2.43) 



92 



The gf at the messenger scale can then be obtained via Eq. (2.40): 



91 = mi/2/lB^- 



(2.44) 



Taken together, the consistency relations, Eq. (2.39) and Eq. (2.42), provide strong con- 



straints that make it highly unlikely, for example, for a generic flavor-blind or GGM spectrum 
to mimic the CMSSM + NUHM. 
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3 Experimental Prospects for the MSSM RGIs 



3.1 Extracting Soft Masses from Observables 

We now turn to the prospects for measuring the RGIs at the LHC. In practice it is necessary to 
convert the observed pole masses into the soft masses entering into the RGIs at the TeV scale. 
Here we comment only briefly on some relevant features of the analysis in the sfermion sector. 
For comprehensive studies, see Ref. [i6] - H8] . 

For the flrst and second squark generations, the off-diagonal components of the mass ma- 
trices are proportional to the corresponding fermion masses and can therefore be neglected. 
Furthermore, the soft masses for the up-type and down-type left-chiral states are equal by 
gauge invariance. This allows one squark pole mass to be written in terms of the others and 
thus removed from the RGIs. For instance, 

m?^ = + cos 2/3 , (3.45) 

where we have ignored the quark masses and the subscript on the squarks refers to the pre- 
dominantly left-handed (and therefore chargino-interacting) mass eigenstate of the squarks in 
the spectrum. 

For the third generation squarks, the off-diagonal components cannot be ignored. Therefore, 
in addition to pole masses, mixing angles must be measured in order to extract the soft masses. 
However, there is still one more observable than there are parameters in each sector, allowing 
the sbottom mixing angle to be removed 

m? cos^ 9t + m? sin^ 9f — m? — m? -|- m? — cos 28 
cos^ Or = ^ , . (3.46) 

mi — m~ 

Here and below, the subscripts 1 and 2 denote the mass eigenstates. This leaves only the stop 
mixing angle to be determined experimentally. 

For the slepton sector, the result depends on the treatment of the sneutrinos. In GGM 
models, the right handed sneutrinos, if present in the spectrum, do not receive SUSY-breaking 
masses at the messenger scale. Their masses are therefore controlled by the Major ana mass 
scale, M/j, of the right-handed neutrino partners, which only couple with the rest of the ob- 
servable sector via the neutrino Yukawa coupling^ Since the neutrino masses 

are smaller or of the order of 1 eV, these Yukawa coupling effects are very small provided 



^In more general flavor-blind models, the SUSY-breaking masses for the right handed sneutrinos are naturally 
of the order of the soft scale, and the same conclusion holds. 
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Mr < O (10^° GeV) . For values of Mr > 10^° GeV, the Yukawa effects must be taken into 
account. However, they will only have an impact on the soft masses if the messenger scale is 
larger than Mr, since otherwise the right-handed sneutrinos decouple at the messenger scale. 
In this work we will make the assumption that either Mr < O (10^° GeV) or Mr > M. This 
justifies our neglect of any right-handed sneutrino effects in the RGIs, as we have implicitly 
done in our treatment of the lepton sector in this article. With this assumption, the mixing 
angles in the sneutrino sector are negligibly small for all generations. 

For the first two generations, as in the squark sector we can remove either a slepton or a 
sneutrino mass. For the third generation we can remove either m?^ or the stau mixing angle 
using the relation 

m?^ = cos^ 9f + m?^ sin^ 9f + cos 2/3 , (3.48) 

where we have ignored the explicit dependence on the relatively small r mass. 

Even if a large number of masses can be determined at the LHC, a finite set of ambiguities 
will arise and affect the calculation of the RGIs. For the first and second generation sfermions, 
unless the mass eigenvalues are nearly degenerate, their measurement is not sufficient to com- 
pute the RGIs, because it must still be determined which eigenstate is the superpartner of 
the left chiral and the two right chiral quarks. An incorrect assignment leads (in the limit 
of vanishing D-terms in the mass matrix), for instance, to a reversal of the soft masses 
and m| for the case of the up-type first and second generation squarks, and these parameters 
do not enter into the RGIs in a symmetric way. Similar reversals may occur for the cases of 
the down- type squarks and the sleptons. Measurements in the chargino sector may be able to 
unravel this ambiguity, since only the left-handed first and second generation sparticles couple 
to the charginos. 

Another problem is the difficulty in distinguishing up-type and down-type first and second 
generation sfermions at hadron colliders. Since the dominant production mechanism is through 
the strong force, the electromagnetic charges may be unmeasurable. However, if our intent 
is to test GGM, it is hkely that > m|^, since this relation holds at the input scale if 
^1 > (which is necessary for a nontachyonic selectron mass in the absence of an FI term) 
and is maintained if Dy > (although this condition is not necessary). Alternatively, one 
can test which, if any, of the two squark mass choices leads to a fulfillment of the necessary 
condition of Db^^s — -D^i ~ 0, and take that as the hypothetical correct assignment when using 
the other RGIs to test the GGM scenario further. The contribution of and to Dy^^ is 
quite different and of opposite sign, and therefore a cancelation of this RGI within errors for 
both choices of the squark masses is highly unlikely. 
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Figure 1: Experimental errors in the low scale determination of the RGIs plotted vs. , assum- 
ing 1% uncertainties in all soft masses and scanning over the high scale parameters as in 



Eq. (3.49). The clear correlation seen in these plots demonstrates that although the invari- 
ants involve many masses squared, the total uncertainties are dominated almost entirely 
by the uncertainty in the squark masses. RGIs with small uncertainties are omitted for 
brevity. 
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Figure 2: Experimental errors in the calculated values of the high scale GGM parameters using the 
RGI reconstruction method. Parameters with small uncertainties are omitted for brevity. 
The spread in the calculated A^. parameters is due mainly to 2-loop effects, and can be 
reduced using the technique discussed in Section 3.3. See Fig. [6j 
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3.2 Measurement Uncertainties in RGIs and GGM Parameters 

It is difficult to estimate liow well the physical quantities of the MSSM can be measured at 
the LHC, even if the sparticles are kinematically accessible. Most mass measurement studies 
in the literature have been performed in mSUGRA/CMSSM benchmark scenarios, under the 
assumption that sparticles will be produced in long cascade decays from strongly produced 
squarks and gluinos (for a detailed study in the context of the LHC and a future linear collider, 
see |100j ). Endpoint methods and other kinematic variables may then be used to extract the 
pole masses. As for the stop mixing angle, in Ref. |101j it was shown that the angle and the 
light stop mass may be inferred from a fit to measurements of ratios of stop branching ratios 
into charginos and neutralinos (for other stop mixing studies, see |1U2[ I1U3] ). The authors 
demonstrated that this method can be effective if the decays are kinematically allowed and if 
the parameters of the neutralino/chargino sector have already been measured. However, it is 
clear that these methods will only be possible in a portion of the allowed parameter space, even 
restricting just to the regions compatible with GGM. Therefore, it is challenging to conjecture 
what can be done at the LHC without reanalyzing the reach point by point in parameter space, 
an exercise which is beyond the scope of this work. 

With this in mind we forgo the application of projected measurement uncertainties to the 
physical parameters and simply apply a universal uncertainty to the soft masses entering into 
the RGIs. This allows us to get a general notion of how well the low scale parameters will have 
to be measured in order to determine the RGIs and the high scale parameters with a certain 
precision. 

We perform a numerical scan over messenger scale inputs, restricting for the purposes of 
illustration to the parameter space of GGM and taking the following ranges of values: 



0.1 < 




< 1.0 TeV^ 


< 




< 1.0 TeV^ 


0.1 < 


MBr 


< 1.0 TeV ; 


2 < 


tan P 


< 50 ; 




M 


< 10^^ GeV 



(3.49) 



For each point in our scanned parameter space of models, we compute the soft spectrum, 
run the soft masses down to the scale fit = 1 TeV, and calculate the resulting central value and 
uncertainties for the RGIs. We then reconstruct the high scale GGM parameters, assuming 
a universal 1% measurement error on the soft parameters. We enforce the conditions that 
the low energy spectrum is within the reach of the LHC and satisfies rriweak > 100 GeV and 
TTT'strong > 300 GcV, whcre the rriweak are the masses of the weakly-interacting sparticles and 
'^strong are the colored sparticle masses. 
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As an approximation, we consider the propagation of errors in quadrature, although this 
is only an estimate since most masses would be obtained in decay chains and therefore their 
experimental errors are expected to be highly correlated. The quadrature combination also 
allows the errors from a larger universal soft mass fractional error to be obtained simply by 
re-scaling the results presented here for 1%. 

Fig. jlj gives the net experimental error in the RGIs as a function of m| . We omit -Dlis 5 
Ib^i and J^^, which have relatively small uncertainties. With the exception of Imzi the 
errors reflect a small y-intercept and small variation for a given m?- . This demonstrates that 
although the invariants involve many masses squared, the total uncertainties are dominated 
almost entirely by the uncertainty in the squark masses, assuming all individual mass errors 
are of the same order. Since in the GGM spectrum the squark mass splittings are proportional 
to smaller gauge couplings, they are typically small compared to the masses themselves unless 
^1,2 ^ ^3- Thus the expected error tends to be controlled essentially by a single parameter, 
as reflected in the scan. In the case of Imz the gluino error must also be included. 

Going a step further, we can estimate the uncertainties in the reconstructed values for 
the high scale GGM parameters. For the this reconstruction relies on the premise that 
5u 7^ within experimental errors, which can be determined by testing the consistency of 
Dyy^jj with zero. If 5u = 5d within errors, then the extraction of gi at the high scale using 
Eq. (2.29) leads to potentially meaningless results. Fig. [2] gives the errors in the calculated 
values of the A.^, 5u,d, and gl at the messenger scale, for those points in the scan satisfying 
I-^>i3hI > '^^^izh' "^^^ brevity we omit the MBr, which have small uncertainties controlled 
entirely by the corresponding gaugino masses. 

It is clear from the range of the y-axes in Fig. [2] that even a precise determination of 
the soft masses can lead to moderate uncertainties in the reconstructed Ar. This is mostly 
due to cancelations between large squark mass parameters, which reduce the value of the Ja/^ 
while increasing the magnitude of the uncertainties. This explains why the value of A3 can 
be determined with relatively good precision, while even for large values of Ai, its uncertainty 
may be as large as its calculated value. In most cases, however, the uncertainties in the Ar are 
smaller than their calculated values, indicating that useful information can be obtained about 
the allowed range of these parameters. The considerable spread in the reconstructed values of 
the Ar is due to 2-loop effects and will be mostly compensated by methods to be discussed in 
the next section. 

As mentioned previously, if measurements suggest 6u = within errors, one can still obtain 
constraints on the Ar by making educated guesses for the gauge couplings at the messenger scale. 
Taking the messenger scale to lie between 10^ and 10^^ GeV, the gauge coupling dependence 
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entering into Eq. (2.32) can then be estimated af 



^ 10 ± 6 , 
^2"^ ^ 4.5 ± 0.5 , 

^3-^^2.5^1.5. (3.50) 

These errors are not very different from what is obtained by propagating soft mass experimental 
errors through the gi reconstruction in the case 5„ 7^ 5d- Observe that the values of the couplings 
are correlated and should satisfy the Ig^ 3 RGI constraints. 

Therefore, one can always use the RGIs to obtain interesting constraints on the A^, even if 
8u = within uncertainties. 

3.3 2-loop Effects in GGM 

One can easily check that the RGIs discussed are not strictly preserved at the 2-loop level and 
therefore it is important to estimate the magnitude of the 2-loop effects. For the scan discussed 
in the previous section, we implemented full 2-loop RG equations for the soft SUSY-breaking 
parameters and the gauge and Yukawa couplings |104j into Mathematica. We compared the 
running of the parameters calculated by our code with those obtained from the public program 
SuSpect |105] . finding excellent agreement. 

The spread in the reconstructed values of the Aj. in Fig. [2] demonstrates that apart from 
experimental error, there is also a theoretical uncertainty in the calculated low scale value 
of invariants since the RGIs defined above have vanishing /3-functions only at 1-loop order. 
The difference between the low and high scale central values of the RGIs and the difference 
between the input and reconstructed central values of the GGM parameters gives the 2-loop 
contributions to these quantities. 

In Appendix A we list the 2-loop /3-functions for all the RGIs, ignoring the small contribution 
from the hypercharge couplings as well as the terms proportional to the trilinear coupling, which 
however have been included in our numerical work. We see that the contributions to Dl-^^ and 
I Ml vanish in the limit that ?/e — ^ 0, and so in general these functions behave as approximate 
RGIs at the 2-loop level. Furthermore, the /3-function of is proportional to the square of 
the bottom Yukawa coupling and it is therefore only relevant for large values of tan /3. Finally, 
the /3-function of Dsy^ does not contain any strong gauge couphng contribution and thus tends 
to be small for low values of the messenger scale. 

To analyze the invariants with larger 2-loop contributions, it is useful to consider a further 
limit of the /3-functions defined in Appendix A in which we turn off the slepton masses, the 

^In principle, if the messenger scale is < 10^ GeV, it may be determined from decays of the NLSP to the 
gravitino inside the detector. Assuming this is the case one could take the messenger scale between 10^ and 
IQi^ GeV, with gj"* « 9 ± 5, g^^ « 4.5 ± 0.5, and g^^ « 2.5 T 1- 
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Figure 3: Left: 2-loop shifts in the low scale values of the RGIs. The x-axis is a simple approximation 
to the 2-loop /3-function. The red line denotes the best-fit for A in terms of our approximate 
/3 function. Right: The ratio of the residual 2-loop effects, after subtracting the best fit /3 
from A, to the experimental errors. We assume 1% uncertainties in all of the soft masses. 
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Figure 4: Left: 2-loop shifts in the low scale values of the RGIs. The x-axis is a simple approximation 
to the 2-loop /3-function. The red line denotes the best-fit for A in terms of our approximate 
/3 function. Right: The ratio of the residual 2-loop effects, after subtracting the best fit /3 
from A, to the experimental errors. We assume 1% uncertainties in all of the soft masses. 
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Figure 5: Left: 2-loop shifts in the low scale values of the RGIs. The x-axis is a simple approximation 
to the 2-loop /3-function. The red line denotes the best-fit for A in terms of our approximate 
/3 function. Right: The ratio of the residual 2-loop effects, after subtracting the best fit f3 
from A, to the experimental errors. We assume 1% uncertainties in all of the soft masses. 
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Figure 6: Experimental errors in the calculated values of the high scale GGM parameters using the 
RGI reconstruction method. Parameters with small uncertainties are omitted for brevity. 
These plots differ from those of Fig. [2] in that the bulk of the 2-loop corrections have been 
accounted for via simple shifts in the RGIs, as detailed in Section 3.3. 
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Figure 7: Percent deviation in left-iianded squark masses necessary to generate values of D (blue 
points) or D^^ (red points) more than la away from zero, assuming 1% experimental 
uncertainties in the soft masses. 
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Figure 8: Percent deviation in right-handed up-type squark masses necessary to generate values 
of -Dbi3 (blue points) or D-^^ (red points) more than la away from zero, assuming 1% 
experimental uncertainties in the soft masses. 
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Figure 9: Percent deviation in right-handed down-type squark masses necessary to generate values 
of Dbi2, (blue points) or D^-^ (red points) more than la away from zero, assuming 1% 
experimental uncertainties in the soft masses. 
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Figure 10: Percent deviation in left-handed slepton masses necessary to generate values of D^-^^ more 
than la away from zero, assuming 1% experimental uncertainties in the soft masses. 
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Figure 11: Percent deviation in right-handed slepton masses necessary to generate values of -Di^g 
more than la away from zero, assuming 1% experimental uncertainties in the soft masses. 



bino mass, and the Higgs masses, set the masses of all squarks equal to mj^, and i/u = a/3/4. 
This reduces the 2-loop /3-functions to simple functions of tan /3, M2, M3, and m^_^. The 2-loop 
corrections can then be estimated from a limited number of parameters as the values of these 
functions multiplied by log(M/TeV). 

To demonstrate the utility of these estimates, we analyze numerically the shifts in the 
RGIs induced by the full 2-loop RGEs. In the left-hand column of Figs. [3]j5j we plot the 
2-loop corrections, A(RGI), against the approximate 2-loop /3-functions. The best-fit line 
passing through the origin then gives an estimate for the average value, (A), of the 2-loop 
corrections. The slope of the line corresponds approximately to log(Mo/TeV), where Mq is 
some intermediate messenger scale. We tabulate the slopes in Table |4} The deviation from 
this approximation is due to residual 2-loop running as well as to the different values of the 
parameters over which we scan. In the right-hand column we subtract off the approximated 
2-loop /3-function shifts, (A)s, from the 2-loop contributions. As, at every point in the scan 
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Table 4: Equations for the best fit line for the 2-loop beta functions, corresponding to Figs, sjjs 



and plot the ratio of A — (A) to the experimental uncertainties against the approximate 2-loop 
/3-functions. The spread in the y-axis is due to the residual 2-loop effects not accounted for by 
the (A)s. 

From the numerical simulation we see that in most cases, even with an optimistic exper- 
imental error of 1% in the soft masses, the experimental error tends to be larger than the 
residual 2-loop effects on the RGIs once the shift is performed to remove the bulk of the 2-loop 
corrections. In those cases it is justified to treat the 1-loop RGIs as true invariants in the de- 
termination of the parameters of the GGM models. For 1^2, the 2-loop corrections can become 
of the same order as the 1% experimental errors. One could then in principle combine the 
uncertainty in due to the 2-loop effects in quadrature with the experimental errors. Notice 
that the determination of the messenger scale by any independent method can serve to further 
reduce most of the uncertainty associated with the 2-loop evolution of the soft parameters. 

In Fig. [6] we repeat the plots of Fig. |2] using the shifted RGIs to compute the GGM input 
parameters. This further demonstrates the advantage of the simple shifts and the subsequent 
dominance of experimental errors over residual 2-loop corrections. 

Finally, we consider more carefully -D^ia, D^-^.^, and -D^i- '^^^ 1^^^^ expected to vanish 

for all flavor-blind models and all three vanish in GGM. To test the power of these invariants 
as discriminants we calculate the percent deviation in the soft SUSY-breaking parameters that 
would lead to a departure from zero. For instance, even a 10% shift in entering into the 
D^-^^ RGI is enough to shift it by more than 5cr outside the range experimentally consistent 
with zero. Similar conclusions hold for Dbi^ and -Dli^. For a few of the soft masses we plot 
in Figs. [7 11 the percent deviation in them that would cause the invariants to take values 
more than la away from zero. Dy^^ and Dbi-^ are quite sensitive to the squark masses; they 
are less sensitive to the generically smaller slepton masses, but for those -Diig is an excellent 
discriminator. As mentioned before, the effect of a larger experimental error in the soft masses 
can be simply seen as a rescaling of the y-axis of Figs. [7 11 
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4 Conclusions 



Low energy supersymmetry provides a predictive framework of physics at the weak scale. How- 
ever, the precise spectrum of the new particles depends on the soft SUSY-breaking parameters, 
which proceed from the unknown mechanism mediating SUSY-breaking from the hidden sector 
to the MSSM at high energies. Knowledge of the soft parameters at the messenger scale will 
greatly contribute to the understanding of UV physics. In principal the RG evolution of the soft 
parameters will allow the determination of their high scale values as a function of the measured 
values at the weak scale. The running of each scalar parameter, however, is complicated by the 
dependence of the ^-function on nearly all other running soft parameters. 

In this work we have proposed the use of Renormalization Group invariants to resolve the soft 
SUSY-breaking parameters at the messenger scale. Of the fourteen RGIs we have considered, 
two may be used to test high scale flavor universahty of the soft parameters, and twelve encode 
information about the thirteen variables associated with the most general CP-conserving flavor- 
universal models. Indeed, if the messenger scale can be determined by alternative methods, the 
whole spectrum at high energies could be established. 

Specific models of SUSY-breaking, however, lead to relationships between the different spar- 
ticle masses at the messenger scale. For example, in General Gauge Mediation, even assuming 
a modification of the soft parameters in the Higgs sector, only eleven free parameters remain. 
Moreover, an additional RGI must vanish and may be used, together with the two RGIs re- 
lated to the fiavor independence, to test consistency of the low scale spectrum with GGM. 
As we have shown in this work, the deviation from zero of these invariants is a very powerful 
discriminant for these models. The remaining eleven RGIs can be used to extract most, and 
in some cases all, of the GGM parameters at the messenger scale. More minimal models, for 
example the CMSSM -|- NUHM, depend on fewer parameters to define the high energy sparticle 
spectrum. In those cases the system is over-constrained, leading not only to the determination 
of all parameters but also to consistency relations between the values of several RGIs. 

Although RG invariance holds only at the 1-loop level, we have shown that in general 
the 2-loop evolution leads to modifications of the RGIs that are smaller than their expected 
experimental errors obtained from aggressive 1% uncertainties in the low scale soft masses. 
Moreover, we have shown that a simple approximation based on the 2-loop ^-functions and 
dependent on a few low-scale parameters effectively describes the 2-loop corrections to the 
invariants. These functions can be directly subtracted from the low energy measured values 
of the RGIs, thereby reducing the theoretical uncertainty. The remaining uncertainty is then 
primarily due to the unknown messenger scale and can be reduced once further constraints are 
obtained on the value of this scale. 

RGIs offer a simple and modular approach to the reconstruction of messenger scale physics. 
The dependence of high scale parameters on observed values at the weak scale is reduced from 
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a set of integrals to a set of algebraic equations. Moreover, each RGl depends only on a subset 
of the low scale masses. Although we have not investigated this possibility in detail, it would 
be interesting to consider strategies for the use of RGIs if only a limited set of the masses 
are extracted at the LHC or a future linear collider. In such a situation traditional methods 
relying on direct integration of the RG equations will fail in the sfermion sector, while if any 
RGIs depend only on the known masses, they may still be used to provide constraints on the 
high scale structure. It would also be interesting to examine in greater detail the use and 
effectiveness of the RGI method in other proposed models for the mediation of SUSY-breaking, 
as well as in non-minimal models of low scale supersymmetry. Additionally, in this work 
we have taken a simple quadrature sum in the estimation of experimental errors, neglecting 
correlations which in many cases may be significant. Finally, a more realistic determination of 
the expected experimental uncertainty in the extraction of the soft breaking parameters could 
alter the overall uncertainties in parameter reconstructions. We leave these investigations to 
future work. 
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APPENDIX A 



2-Loop /3-Functions for the 1-Loop RG Invariants 

In the following set of 2-loop ^-functions, we set the hypercharge gauge coupling gi and the 
soft trilinear couplings Ai to zero. The subscripts on the (3 functions are in correspondence 
with the invariants defined in the text. 



= ^ (3^2 (rul + + ml^ + 2M|^ + 2 (^m| + + m^^^ 

+ (3^72 (rnl^ + ml + + 2Ml^ - 2 (m\ + 2m|^ + ml^ + m\ + m|„^ y^\j yl 
+ 2 {ml^ + + m|„^ |/^[^ , (A.l) 
^ (3^72' (m^3 + + m^^ + 2M2^) yl + 2 (^m?3 + + m^^) |/,^) , (A.2) 
(3^?2^?3 (^fl + M2M3 + ^ , (A.3) 
^D,,3, = ^ (3^2 (m| + M2M3 + Ml^ - 5 (^m| + + mf,, + 2M|^ yl 

- 5 (^m?3 + + + 2M|^ yl + 10 (^m|^ + ml^ + m^^ + 2M2^ yl^ 

+ 5 (^8m|^^32y2 + ^lm]jjjl + ley^M^yJ ^ 4^|^^4 ^ 4^2^^^4 ^ 3^2^^2^2 ^ 3^2^^2^2 



^0^1 - 



'^3 (^^^3 + ml + + 2M|^ + (^2mJ^ + m|3 + m\ + m^^^ yl 

- 2 (^m|^ + m?3 + m^^^ + m\yl {^gl + 4?/^ + 'iyl - ) ) , (A.4; 

^l/d (9^2^^ - ^^9lm\^ + 18<72'M2^ - ZlglMl - ^m\f, - 6 [m\ + + 2mi 

+ m\ {^gl - 2 {%gl + 3 (^y^ + y,') ) ) + ^^3 (^Q^^' - 2 (^8^| + 3 (^y^ + y,^) ) ) 

fA.5' 
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6407r4 
1 

128^ 
1 



+ 2 + m?3 + ml}j yt + m\yl - - - ^yl + t/^^ 

+ [yl - 2yl}j + 2 (^gl (^Ml + M2M3 + + M^^ (^yj + ^2 _ 3^2^ ^ 

(A.6) 

27glM2 + W^M^ , (A.7) 



25^2^2 + 2AglAhj , (A.8) 
g^^^Ms + UglM^ , (A.9) 



5/ 



, (A.IO) 
]^^^2 (^"^!i^2 + ^l^gl + ^l^al + ^rn^Q^d + ^TT^l^d + 3^13^2 + ^2^L + dWn^ 
+ lllglMl + A9>glMi + ASgjM^Ms - UM^yl - AM^yl - 12M|y^^ , (A.ll) 
Y^^3 (3m|y| + 3m\gl + 3m|^| + 6m|^^| + 6m|^^| + 6m|^^2 ^ 3^2^ ^2 

+ Sml^^l + Smlgl - SOglM^ - ISglM^ (^M^ + Mg^ + SM^yl + SM^yl^ , 

(A.12) 

1 



3207r4 
1 

320^ 



399^2 + 352^2 - 92yl - 24yl - 86ylj , (A. 13) 

- 63^2' - 121^1 + 29yl + 9yl + Sbyl] . (A.14) 
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